Control of gradient-driven instabilities using shear Alfven beat waves 
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A new technique for manipulation and control of gradient-driven instabilities through nonlinear 
interaction with Alfven waves in a laboratory plasma is presented. A narrow, field-aligned density 
depletion is created in the Large Plasma Device (LAPD), resulting in coherent, unstable fluctuations 
on the periphery of the depletion. Two independent shear Alfven waves are launched along the 
depletion at separate frequencies, creating a nonlinear beat- wave response at or near the frequency 
of the original instability. When the beat- wave has sufficient amplitude, the original unstable mode 
is suppressed, leaving only the beat- wave response, generally at lower amplitude. 

PACS numbers: 52.35.Bj, 52.35.Mw, 52.35. Qz 



In magnetized plasmas in the laboratory, the loss of 
heat, particles, and momentum across the confining mag- 
netic field is predominantly caused by turbulent trans- 
port associated with pressure-gradient-driven instabili- 
ties [U [2] . Controlling these instabilities to mitigate the 
transport they produce is therefore highly desirable for 
magnetic confinement fusion devices such as tokamaks. 
Direct control of gradient-driven instabilities has been 
achieved using arrays of electrodes or external saddle coils 
to induce fluctuating parallel currents in a cylindrical 
plasma device [3, 4 . In these experiments, when the ex- 
cited electric fields matched the frequency and wavenum- 
ber of a drift-wave mode in the device the observed 
turbulent spectrum collapsed onto the single coherent 
driven mode (synchronization) and turbulent transport 
is reduced [5 . Extending the transport control enabled 
through mode-synchronization to high-temperature fu- 
sion devices is highly desirable but it would not be possi- 
ble using material electrodes. It is also not clear whether 
saddle coils would be effective at driving currents for syn- 
chronization in a high-performance device. 

Turbulence control with core plasma access might be 
achieved using externally-launched radiofrequency (RF) 
waves. High-power RF is expected to be an important 
tool for heating and current drive on next-step fusion ex- 
periments such as ITER [6 and could be available for 
application to turbulence and transport control. Cre- 
ation of flow shear for turbulence suppression [7 using 
RF waves has been theorized [8 and recently demon- 
strated in the Alcator C-Mod tokamak [9l [10]. Direct 
modification of gradient-driven instabilities is also possi- 
ble through interaction with an RF field. Control of drift- 
wave fluctuations was demonstrated in Q-machine exper- 
iments using lower hybrid waves to affect plasma proper- 
ties and the drift- wave dispersion to suppress growth [11 - 
[T4] . The modification of other transport related instabil- 
ities, such as the ion temperature gradient (ITG) insta- 
bility, in the presence of RF driven fast magnetosonic 
waves has also been suggested and investigated theoreti- 
cally [ig. 

This Letter reports on novel laboratory experiments 



using RF waves, in particular shear Alfven waves 
(SAWs), to modify gradient-driven instabilities. Two co- 
propagating SAWs with slightly differing frequency are 
launched along a narrow field-aligned density depletion 
on which coherent (dominant single-frequency) unstable 
fluctuations are observed. Modification and control of 
the instability is seen when the nonlinear beat response 
exceeds an amplitude threshold and the frequency of the 
beat wave is near the frequency of the gradient-driven 
instability. In this case, the original unstable mode is 
reduced in amplitude or suppressed completely, leaving 
only the beat-driven response, generally at lower am- 
plitude. A reduction in broadband fluctuations is also 
observed accompanying the suppression of the primary 
mode. This control is observed for a finite range of beat 
wave frequencies around the unstable mode frequency; 
this range widens with increased beat wave amplitude. 

This research was carried out on the upgraded Large 
Plasma Device (LAPD) at the Basic Plasma Science 
Facility (BAPSF) at UCLA jjG,. The LAPD creates 
a 17.5m long, '^0.6m diameter cylindrical magnetized 
plasma. The plasma is created at 1 Hz using an 
indirectly-heated barium oxide coated cathode discharge 
(discharge length is typically 10-20ms). Bulk parameters 
for the helium plasmas used in these experiments are: 
He - lO^^cm-^, Te - 6eV, and B = 800G. Ion tempera- 
ture in LAPD is typically Ti ^ leV [18 , although it was 
not measured directly in these experiments. Under these 
conditions, the electron thermal speed is greater than the 
Alfven wave speed (the Kinetic Alfven Wave regime p!7]). 

In these experiments, temperature, density and float- 
ing potential are determined using a triple Langmuir 
probe. The temperature and density measurements have 
been quantitatively verified using microwave interferom- 
etry and swept Langmuir probe measurements. Mag- 
netic field fluctuations were measured using magnetic 
pickup coils, differentially wound to eliminate electro- 
static pickup. The probes are mounted on motorized, 
computer controlled probe drives. Two-dimensional aver- 
age profiles are acquired over many repeatable discharges 
by moving a single probe shot-to-shot. Two-dimensional 



cross-correlation functions are achieved using a second, 
stationary reference probe separated axially (along the 
background field) from the moving probe. 

Conditions are established for the instability formation 
through creation of a field-aligned axial density depletion 
embedded in the center of the LAPD plasma column. 
This type of condition has been extensively investigated 
in the LAPD, and the presence of gradient-driven drift- 
Alfven-wave instabilities in such conditions is well docu- 
mented [191 [in] • The depletion is formed by inserting a 
small metal blocking disk between the anode and cath- 
ode, selectively blocking the high energy primary elec- 
trons that ionize the neutral gas and leaving a lower tem- 
perature and density striation along the central axis of 
the machine. As shown in Figure Ilia) the pressure drops 
by about 50% in the depletion, which is ~6 cm wide. 
The measurements shown are of ion saturation current 
(/sat oc ne\/Te) and pressure Pg, with Pg calculated from 
Tg and /gat- Figure [lib) shows measured /gat fiuctua- 
tions, which are localized to the density gradient associ- 
ated with the depletion. Figure md) shows the spatially- 
averaged fast Fourier transform (FFT) power spectrum 
of /sat, indicating a narrow-band spectrum with power at 
~ 6 kHz and harmonics. Cross-correlation measurements 
shown in Figurelllc) of the spontaneous fiuctuations indi- 
cate a m = — 1 azimuthal mode number for the dominant 
mode, with propagation in the electron diamagnetic di- 
rection. Many of the measured properties of the unstable 
modes (frequency, direction and speed of propagation) 
agree with predictions for drift waves. However, Fig- 
ure file) shows that the measured phase between density 
fiuctuations and potential fiuctuations {Ox) is approxi- 
mately 180°, which is inconsistent with the {)° < 6^ < 90° 
predictions and observations of saturated, single mode 
drift waves [21 . 

It is important to note that in these experimental con- 
ditions there are two sources of free energy, the pres- 
sure gradient and an imposed potential gradient from 
the blocking disk, resulting in two separate fiows from 
the E X B drift and the V x P drift. Thus there are 
two possible instabilities, the drift-wave instability, re- 
sulting in drift- waves with axial mode number ny > 0.5, 
and the Kelvin-Helmholtz instability, resulting in "fiute- 
like" modes with ny =0. Numerical eigenmode solutions 
to the linearized Braginskii two-fiuid equations [22] were 
performed, using experimentally measured density, tem- 
perature, and potential profiles. These solutions indicate 
that both types of modes are linearly unstable and have 
real frequency approximately equal to the measured fre- 
quency, leaving open the possibility that the saturated 
single unstable mode observed might be of either type of 
instability or a hybrid mode driven by both sources of 
free energy. 

A two-strap "picture-frame" style antenna [23l [24] is 
inserted into the plasma at the far end of the deple- 
tion, 13 m from the plasma source. The antenna is 




FIG. 1. [Color] (a) Contours of measured plasma pres- 
sure in the depletion. (b) Contours of /sat fluctuation 
power, with pressure contours overlaid in white, (c) Cross- 
correlation showing m=-l structure of mode, (d) Spatially 
averaged hat FFT spectrum, (e) Spatially averaged cross- 
phase (^x) between /sat and Vf fluctuations 



centered on the depletion and aligned to launch SAWs 
along the depletion. Each strap of the antenna is driven 
and controlled independently, giving rise to two perpen- 
dicularly polarized, co-propagating SAWs. The frequen- 
cies of the two waves are offset by a few kHz, and with 
each frequency near half of the ion cyclotron frequency 
(/ci = 304 kHz at 800G). The peak wave magnetic field 
for SAWs used in this study ranges from O.IG to IG 
(1 X 10-"^ < SB/B < 1 X 10-^). In the absence of the den- 
sity depletion, the two SAWs interact nonlinear ly, driving 
a nonresonant quasimode at their beat frequency [25 . 

To investigate the interaction between the beating 
SAWs and the instability driven by the density deple- 
tion, a series of experiments were performed in which 
the beat frequency was set near the instability frequency 
and varied shot-to-shot. Fig. |2la) shows the measured 
^sat power spectrum near the depletion gradient maxi- 
mum (at approximately x=2 cm, y=-l cm (see Fig. fl] 
axes) 6 m from the antenna and 7 meters from cathode) 
as the beat frequency is varied over < /beat < 15 kHz. 

The beat-driven response is observed to be strongest 
when the beat frequency is close to the instability fre- 
quency (6 kHz), with peaks centered at 5.5 kHz and 
8.25 kHz as shown in Fig. [2lb). In this particular 
dataset, the peak RMS amplitude of the beat-response 
is SIsat/Isat = 15%. The measured beat-wave response is 
significantly stronger than that observed in the absence 
of the depletion, where for the same launched SAW pa- 
rameters, a beat- wave amplitude of ^/sat/^sat ^ 5% is 



observed. 

The most striking feature of Fig. [2] is the effect that 
the SAW beat wave has on the unstable mode at / ~ 6 
kHz. As the beat- wave frequency approaches 6 kHz from 
below, the spontaneous mode begins to decrease in power 
and is observed to downshift in frequency, meeting the 
beat wave where it appears to resonate near 5.5 kHz. As 
the beat wave is driven to higher frequency, the sponta- 
neous mode appears completely absent, as if the growth 
of the mode is suppressed in the presence of the driven re- 
sponse. Fig.[2|b) illustrates the shift of fluctuation power 
from the instability to the driven fluctuations by show- 
ing the power at the unstable mode frequency (6 ± 0.4 
kHz) and power in the beat frequency band {/beat i 0-4 
kHz) as a function of driven beat frequency, essentially a 
horizontal and diagonal line cut of the data in Fig. |2|a). 
Also of note is the apparent reduction in broadband fluc- 
tuations while the unstable mode is suppressed. This is 
made clearer in Fig. [2Vc), showing the measured power 
spectrum in the absence of the SAW beat wave (both 
antennas set to the same frequency, black line), and the 
spectrum in the presence of a 6 kHz (blue line) and an 
8.25 kHz (red line) driven beat wave. As the beat fre- 
quency is increased further, eventually the instability re- 
turns as the beat-driven response appears to diminish. 

The peaks in Fig. [2|b) suggest that the SAW beat 
wave is resonantly driving a linear mode at / ~ 5 kHz 
and / ~ 8 kHz and this conjecture is supported by the 
spatial structure of the beat-wave response found using 
cross-correlation measurements. Figure |3] shows the two- 
dimensional cross-phase (row a) and spatial distribution 
of fluctuation power (row b) for three cases: the spon- 
taneous mode in the absence of SAWs (column 1), the 
beat-wave response when driven at 6 kHz (column 2), 
and the beat wave response when driven at 8 kHz (col- 
umn 3). Note that for the two beat-wave cases, the power 
in the beat-wave response is concentrated on the density 
gradient (Figs. [3]^b2) and (b3)), similar to the sponta- 
neous mode (Fig. IsFbl)). In contrast to this observation, 
in the absence of the depletion, the beat-wave has a much 
broader spatial distribution consistent with the pattern 
of the launched SAWs (not shown here). The correlation 
function of the beat driven response at 6 kHz (Fig. 3 ^a2)) 
is nearly identical to the spontaneous mode (Fig. 3'al)) 
and indicates an tti = —1 mode structure. The beat 
driven response at 6kHz also rotates in the electron dia- 
magnetic drift direction. The peak in the beat response 
at 8kHz corresponds to an ?7i = — 2 mode, as shown in 
Fig. Has). 

The appearance of the control of the instability shown 
in Fig. [2] depends strongly on the strength of the beat- 
wave drive, a parameter controlled by increasing the am- 
plitude of the launched SAWs. Fig. [IJa) shows the vari- 
ation in the beat-wave and the instability amplitudes as 
the SAW antenna current is increased (for fixed 8kHz 
beat frequency). A threshold is clear for affecting the in- 
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FIG. 2. [Color] (Units normalized to power equivalent of max- 
imum RMS fluctuations of 15%) (a) Measured /sat fluctuation 
power as a function of the SAW beat frequency (x axis) and 
FFT frequency bin (y axis), (b) Power at the spontaneous 
mode frequency (/ = 6 ± 0.4 kHz) and at the SAW beat 
frequency (/ = fbeat ± 0.4 kHz), (c) /sat fluctuation power 
spectra at 0.0, 6.0, and 8.25 kHz SAW beat frequency (marked 
with vertical lines on the 2-D bi-spectral plot above in part 



stability, requiring the beat- wave amplitude to be on the 
order of 10% of the amplitude of the original unstable 
mode. It is also interesting to note that the beat-wave 
response saturates as the instability is completely sup- 
pressed. Also shown is the instability power as a function 
of antenna current in the case of zero beat frequency. A 
reduction of the instability power at high SAW amplitude 
is observed which can be attributed to profile modifica- 
tion due to electron heating by the launched SAWs [24 • 
This effect contributes to modification of the instability 
in the presence of SAWs, but is less significant than the 
reduction observed with non-zero beat frequency. 

The width in beat frequency over which control of the 
instability is observed also varies with amplitude of the 
driven beat wave response. Fig. mb) shows the FFT 
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FIG. 3. [Color] 2D phase (from cross-correlation) and spa- 
tial power distribution of /sat fluctuations for the spontaneous 
mode (no SAWs) and 6 kHz and 8 kHz beat waves. 



power at the unstable mode frequency as a function of 
beat wave frequency and antenna current (which controls 
the amplitude of the interacting Alfven waves) . At higher 
currents, the instability (m = —1) begins to be sup- 
pressed in favor of the higher frequency mode driven by 
the beat-wave {m = —2). The white contour in Fig.[4Fb) 
indicates where the instability power is at 50% of the 
'background' power with the antennas driving a zero beat 
frequency Alfven wave. The broadening of the range over 
which the instability is controlled as the beat-wave power 
is increased is indicated by this contour. Similar obser- 
vations are made in VINETA [26], where the effect was 
compared to "Arnold Tongues" seen in nonlinearly driven 
unstable oscillators [271 EH] • 

The exact mechanism by which the beating SAWs cou- 
ple nonlinearly to and control the unstable mode is still 
under investigation. Experiments have established that 
the parallel wavelength of the beat wave is essential. Beat 
waves generated by co-propagating SAWs have very long 
parallel wavelength, as long or longer than the length 
of LAPD at low frequency [25], similar to expectations 
for the parallel wave number of instabilities in the de- 
pletion. Experiments were performed in which counter- 
propagating SAWs interacted along the depletion, driv- 
ing beat waves with much shorter parallel wavelength. In 
this case, much smaller response at the beat frequency is 
observed and the control behavior shown in Fig. [4] is not 
achieved. 

An important question is whether or not transport is 
affected through this beat- wave control process. The 
unmodified instability has density-potential crossphase 
(Figure [l](e)) that is inconsistent with fluctuation-driven 
particle transport. Past experiments on filamentary pres- 
sure striations with dominant coherent modes in LAPD 
have revealed classical transport levels [29]. In the pres- 
ence of the beat-wave in the experiments reported here, 
the density-potential crossphase is unaltered. Future ex- 
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FIG. 4. (a) Power in the instability band(6.5 =b 0.4 kHz, solid 
line) and driven beat response band (8.0 ± 0.4 kHz, dashed 
line) versus SAW drive current for the case of 8.0 kHz beat 
frequency, along with a reference instability power (6.5 di 0.4 
kHz, dotted line) at zero beat frequency (both antennas at 
same frequency), (b) Fluctuation power in the instability 
band (6.5 di 0.4 kHz) vs beat-frequency vs antenna current. 
The white contour is the 50% suppressed level at each power. 



periments will investigate the interaction between SAW 
beat waves and broadband turbulence in LAPD, to see 
whether or not transport control can be affected, as has 
been seen in mode synchronization experiments on the 
VINETA device [5. 

In summary, a novel approach to control of gradient- 
driven instabilities in magnetized laboratory plasmas has 
been presented. Beating SAWs have been observed to 
nonlinearly interact with and control unstable fluctua- 
tions on a filamentary density depletion in LAPD. These 
results provide motivation for studying the interaction 
between RF and instabilities in other magnetized plas- 
mas, including fusion plasmas, e.g. ICRF beat-wave in- 
teraction with microturbulence. ICRF beat waves have 
been successfully used to produce measurable plasma re- 
sponse at the beat frequency and to excite Alfven eigen- 
modes in JET and ASDEX Upgrade [30l|3l], motivating 
future experiments in which the beat wave frequency is 
lowered into the range of drift-wave turbulence. 
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